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Abstract 

An Eulerian-Lagrangian approach to incompressible fluids that is con- 
venient for both analysis and physics is presented. Bounds on burning 
rates in combustion and heat transfer in convection are discussed, as well 
as results concerning spectra. 

Incompressible fluids are described by the Navier-Stokes equation. Turbu- 
lence (0, 1^, jsj) experiments provide measurements that correspond to aver- 
ages of certain quantities associated to the variables appearing in the Navier- 
Stokes equation. The present mathematical knowledge about the Navier-Stokes 
equations is incomplete. Some of the quantities measured in experiments are 
accessible to mathematical theory. They are usually low order, one-point bulk 
averages like the time average of integrals of squares of gradients. Most other 
measured quantities are not amenable to rigorous quantitative a priori analysis. 
Turbulence is concerned with statistical or collective properties of fluids. Never- 
theless, the main impediment to progress in the rigorous analysis of turbulence 
is the present lack of understanding of possible blow up in individual solutions 
of the Euler and Navier-Stokes systems. 

I will discuss briefly the blow up problem and present an Eulerian-Lagrangian 
approach to fluids. I will also give examples of low order one-point bulk quantities 
that can be treated with present knowledge and discuss results on certain two- 
point quantities. 

1 An Eulerian-Lagrangian Approach to Fluids 

I will start by recalling that the Navier-Stokes-Euler system can be written as 
an evolution equation for the three-component velocity vector u = u{x,t), 

du 

— + It • Vm + Vp = uAu + /; 
ot 

the pressure p = p{x, t) preserves incompressibility 

V • M = 0. 



The Euler system is obtained if the kinematic viscosity vanishes, — 0; the 
Navier-Stokes system ii v > 0. Boundary conditions are different for the two 
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systems. The blow up question can be stated in its simplest form for the pure 
initial value problem: are there any smooth initial data with finite energy that 
lead to solutions that diverge in finite time? The answer is not known. 

For a blow up in the Navier-Stokes system one would have to have a finite 
time divergence of an eddy-viscosity-like quantity: 

/ sup \u{x + r,t) ~ u{x, t)\^dt < oo ^ u e C°° . 

Jo x,r 

By contrast, it is known (Q) that 

i-T 

sup \u(x + r, t) — u(x,t)\dt < oo. 



In a situation in which all velocities are finite no singularities can appear in 
the Navier-Stokes equations. One could accept the finiteness of velocities as a 
physical hypothesis. For the Euler system this hypothesis would not be sufficient 
to ensure smoothness of solutions ( |^ ) . A sufficient condition for regularity in the 
Euler equations is the finiteness of the time integral of the maximum magnitude 
of vorticity (|^). The vorticity (the curl of velocity or anti-symmetric part of 
the velocity gradient), w = V x u obeys a quadratic evolution equation. The 
magnitude of vorticity obeys 

Dt\uj\ = a\uj\ 

with Dt — dt + u ■ W the material derivative along particle paths. The scalar 
stretching term a is related to the magnitude of vorticity by a principal value 
singular integral (0): 

a{x)= P.V.J D{y,ax),ax + y))Hx + y)\^ 

\uj{x)\ 

The smooth, mean zero function of three unit vectors D vanishes when two of 
its arguments are on the same line. Consequently, if the direction of vorticity ^ 
is Lipschitz then the singular integral representing a is mild and the solutions 
remain smooth (||], [0). This is a generalization of the two dimensional 
situation where = and the solutions remain smooth. The possibility of 
blow up due to strong vortex stretching is not removed by the previous result; 
the result only precludes blow up in a smooth vortex line field. Blow up can 
occur also because of strain intensification (the strain matrix is the symmetric 
part of the gradient of vorticity). I will present now a description of the Euler 
equations that is convenient for analysis (||ll|) and allows for a clearer geometric 
picture of the possible singularity formation. I will present the results only in 
the periodic case for simplicity of exposition. The case of decay at infinity is 
almost identical. 
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Theorem 1 A function u{x,t) solves the incompressible Euler equations if and 
only if it can be represented in the form u = ua, 

u),i.,t)=u^iAi.,t)f-^:^-'-^ 

where A{x,t) solves the active vector equation 

dA{x,t) 
di 

with initial data 



+ UA{x,t) ■ S/A{x,t) = 



A{x,0) = X. 

uo is the initial velocity and nA{x,t) is determined up to additive constants by 
the requirem.ent of incompressibility, V ■ ua = 0. 

Let us denote by 

P = (1- VA^V-) 

the Leray-Hodge projector on divergence-free vectors. The local existence theo- 
rem requires just one derivative to be Holder continuous: 

Theorem 2 Let uo be a divergence free C^ '^ periodic vector valued function of 
three variables. There exists a time interval [0, T] and a unique C([0, T]; C^'"^) 
spatially periodic vector valued function S{x,t) such that 

A{x,t) = x + S{x,t) 

solves the active vector formulation of the Euler equations, 

dA 

— +UA-VA = 0, 

ot 

UA^^{{VA{x,t)yuo{A{x,t))} 

with initial datum A{x, 0) = x. 

The proof of this result is based on an identity that removes the apparent 
ill-posedness, on singular integral calculus, and on the use of the method of 
characteristics. In the active vector formulation, the " back-to-labcls" map A 
has conserved distribution. Its time evolution is a smooth, volume-preserving 
rearrangement. Singularities can occur only if the gradient map VA diverges 
rapidly in finite time: 



|2, 



/ sup|Vy4(x,t)|^dt<oo^ C°°. 



Thus, would-be singularities arc gradient singularities in a conserved quan- 
tity, similar to shocks in conservation laws, but with the significant difference 
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that the characteristic flow is measure preserving. The formula relating the ve- 
locity to the value of the spatial gradient at the same instance of time (p2[, 
[ p^ , [ p^ , p^) may have important mathematical consequences. On one hand, 
conservation of kinetic energy confers a constraint to the growth of VA. On 
the other hand, the formula suggests that near regions of high gradient of A the 
velocity is exceedingly high, making the regions of high gradient difficult to track 
and perhaps unstable. Dynamical stability or instability of blow up modalities is 
a difficult subject. There are obvious space and time symmetries (for instance, 
a minute delay of blow up), that clearly should not be categorized as instabil- 
ities. Nevertheless, even relatively simple PDEs can exhibit the coexistence of 
a variety of dynamical behaviors, including several stable blow up modalities, 
stable time independent solutions, unstable blow up modalities and dynamical 
connections between the unstable behaviors and the stable ones ([@). 

I will pass now from the blow up problem to some more tractable questions 
about average properties. One can obtain rigorous upper bounds for certain 
bulk averages of solutions of Navier-Stokes equations. Lower bounds are harder 
to obtain. Upper bounds for bulk averages of low order moments for Rayleigh- 
Benard convection will be described further below; I start with a lower bound 
for the burning rate in a simple model of turbulent combustion. 

2 Bulk Burning Rate 

Mixtures of reactants may interact in a burning region that has a rather com- 
plicated spatial structure but is thin across (iQ)- This reaction region moves 
towards the unburned reactants leaving behind the burned ones. When the re- 
actants are carried by an ambient fluid then the burning rate is enhanced. The 
physical reason for the observed speed-up is believed to be that fluid advection 
tends to increase the area available for reaction. What characteristics of the 
ambient fluid flow are responsible for burning rate enhancement? The question 
needs first to be made precise, because the reaction region may be complicated 
and, in general, may move with an ill-defined velocity. An unambiguous quantity 



V representing the bulk burning rate is defined in (|18 ) and explicit estimates 
of V in terms of the magnitude of the advecting velocity and the geometry of 
streamlines are derived. In situations where traveling waves are known to exist, 

V coincides with the traveling wave speed and the estimates thus provide auto- 
matically bounds for the speed of the traveling waves. The main result of ([Q) 
is the identification of a class of flows that are particularly effective in speed- 
ing up the bulk burning rate. The main feature of these "percolating flows" is 
the presence of tubes of streamlines connecting distant regions of burned and 
unburned material. For such flows we obtained an optimal linear enhancement 
bound V > KU where U represents the magnitude of the advecting velocity and 
K is a, proportionality factor that depends on the geometry of streamlines but 
not the speed of the flow. Other flows and in particular cellular flows, which 
have closed streamlines, on the other hand, may produce a weaker enhancement. 



Navier-Stokes Equations and Fluid Turbulence 



5 



The instantaneous bulk burning rate is defined by the formula 

f dT 

^(*)=y -g^{x,y,t)dxdy 

where the integral extends over the spatial domain D, taken here for simplicity 
of exposition to be a two-dimensional strip of unit width and infinite length 

< 2/ < 1, — oo < a; < oo. 

The temperature T is assumed to obey Neumann boundary conditions at the 
finite boundaries and to obey 

T(-oo,t/) = l, T(oo,t/) = 0. 

The simplified model is a passive reactive scalar with a KPP nonlinearity 

Tt + u-VT- kAT = ^T{1 - T). 
with prescribed velocity u that satisfies 



L 



u{x,y,t) dy = 0, V • u = 0. 



The constant vq represents the speed of a stable one-dimensional laminar {u = 0) 
traveling wave. We start with a very general lower bound: 

Theorem 3 For arbitrary initial data obeying 

0<To{x,y) < 1. 
one has the general lower bound 

V{t) > Cvo (l - e"^^ . 

The proof is based on the product lemma 

Lemma 2.1 There exists a constant C > such that 

0<T{x,y) < 1, 

T(-oo,y) = 1, T(oo,t/) = for any y € [0, 1]. 

implies 

Y T(l - T) dxdy^ (^j |VT|2 dxd'^ > C. 
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Although information about the velocity is not present in the general result, it 
nevertheless shows that this model does not permit quenching. Also, the general 
lower bound applies to the homogenized version of the equations as well. 

For a very general class of velocities u{x, y, t) it can be shown that the bulk 
burning rate may not exceed a linear bound in the amplitude of the advecting 
velocity. For a large class of flows we proved lower bounds on the bulk burning 
rate that are linear in the magnitude of advection. We denote by 



the time average of the instantaneous bulk burning rate. The main result of 
([l^) is too technical to state here precisely but its meaning is that presence of 
coherent tubes of streamlines connecting unburned and burned regions enhances 
the burning rate 



as long as the velocity spatial scales are not too small compared to the reaction 
length scale , and the time scale of change of the advecting velocity is not too 

small compared to tq = max[;j^, ^] where H is associated to the width of the 

coherent tubes of streamlines. For instance, a result concerning mean zero shear 
flow of the form 



can be stated as 

Theorem 4 Let us consider an arbitrary partition of the interval [0, 1] into 
subintervals Ij = [cj — hj,Cj + hj] on which u{y) does not change sign. Denote 
by D+ the unions of intervals Ij where u{y) > and u{y) < respectively. 
Then there exist constants C± > 0, independent of the partition, u(i/), and the 
initial dataTQ{x,y), so that the average burning rate {V)t satisfies the following 
estimate: 



T 








{V)r > KU 
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for any r > tq = max ^ . (H = Ij. Here I = k/uq. The constants c± are 
defined by 




The main result of applies to a large class of flows that are not neces- 

sarily spatially periodic, nor shears, and can have completely arbitrary features 
outside the tubes of streamlines. The bulk burning rate is still linear in the 
magnitude of the advecting velocity, no matter what kind of behavior (closed 
streamlines, areas of still fluid, etc.) the flow has outside the tubes. The pro- 
portionality coefficient depends on the geometry of the flow in a rather complex 
manner. These bounds can be extended to larger classes of chemistries. The 
lower bounds, however, are not yet available for models in which there is a feed- 
back coupling of temperature on the velocity. For such models upper bounds can 
be derived. In the next section we will discuss upper bounds for heat transfer. We 
will concentrate on the simplest coupling, mediated by gravity in a Boussinesq 
approximation and discuss the canonical case of Rayleigh-Benard convection. 

3 Bulk Heat Transfer 

The equations for Rayleigh-Benard convection in the Boussinesq approximation 
are 

du 

— — h u ■ Vii + \7p = aAu + aRaeT, 
at 

V • M = 

dT 

+ U-VT ^ AT. 

ot 

The vertical direction e of gravity is singled out. We consider as spatial domain a 
box of height 1 and lateral side L. The velocity vanishes at the boundary, T — 1 
at the bottom boundary and T = at top. The Nusselt number is the space-time 
average of the flux of temperature across horizontal cross-section planes. From 
the equations of motion it follows that 

(|vr|2) = N. 

and also 

(|Vup) = Ra{N ^l). 
Here < • • • > is global space-time average. The general rigorous result here is 

(HI, 0) 

Theorem 5 There exists an absolute constant C , independent of Rayleigh num- 
ber Ra, aspect ratio L and Prandtl number a such that 

N < i + cVml 



holds for all solutions of the Boussinesq equations. 
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The experimental data (|^) point however more to results of the type N ^ 
Ra^ or N ^ Rqt . The exponent i can be obtained rigorously for a simplified 
model. Consider the infinite Prandtl number equations for rotating convection 

{dt+u-V)T = AT 

-Am - E^'^v+px = 
-Av + E^^u + py = 
-Aw+p:, = RT. 

V • u = 

with boundary conditions: p, T) periodic in x and y with period L; 

M, V, and w vanish for z = 0, l,r = Oatz = l,r=latz = 0. One can prove 

(H) 

Theorem 6 There exist absolute constants ci, C4 so that the Nusselt number 
for rotating infinite Prandtl-number convection is bounded by 

N-l< 

min|cii?t; {c2E^ + C3E)R'^;c4R^/^{E-^ + \og^ R)iY 
This coincides, in the limit of no rotation E oo with a logarithmic correc- 



tion (1 23 1) to the I exponent. The bound also shows that strong rotation iJ — > 
stabilizes the system and that increasing rotation may result in a non-monotonic 
behavior of the Nusselt number, as observed in experiments. 

Consider now the horizontal average T(z, t) of T{x, y, z, t) and define 

n = (|v(r-T)p). 

Note that 

n< N 



Theorem 7 For the full Boussinesq system 

N <1 + c{nRa)i 
holds. For the infinite Prandtl number system 

N <1 + C (i?a(log+ Raf^/^) ^ 

holds. 

Note that, if n < is used then the result recovers the exponents ^ for gen- 
eral Rayleigh - Benard and logarithmically corrected i for the infinite Prandtl 
number case. But the rigorous appearance of the exponent | is perhaps not 
coincidental. 
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One of the technical ingredients for the proof of these results concerns zero 
order operators that are not translation invariant 

where w = (A|,^)~-^/ is the solution of 

A^w = f 

with horizontally periodic and vertically Dirichlet and Neumann boundary con- 
ditions w = w' = 0. 

Theorem 8 For any a € (0, 1) there exists a positive constant C'a such that 
every Holder continuous function 6 that is horizontally periodic and vanishes at 
the vertical boundaries satisfies 

llB^IUoo <C„||^|Uoo (l+log+ \\9\\co..)\ 

The proof of this result is based on a pointwise bound on exponential-oscillatory 
sums of the type: 

K{x,z) = e^*'=-^TO^e-^'"'= 

where e = e{z) > 0, ruk = ^\k\ and e{z) = => z = zq- The sum is singular at 
z = zq and the pointwise bound 

p+2 

\K{x,z)\<Cp[\xf + e\z)] ^ 

is obtained using the Poisson summation formula. 

The bounds on bulk one-point quantities presented above are among the 

most successful areas of mathematical and experimental agreement. The reason 
is perhaps that the quantities involved are numbers, albeit numbers depending 
on a parameter. The next step beyond the description of bulk one-point averages 
is the description of power spectra. These are asserted to have some universal 
features in physical turbulence theory; we present some mathematical results in 
the next section. 

4 Spectra 

Unlike bulk one-point quantities, spectra are averages of functions. There are 
some well-established spectra in the physical literature associated to small scale 
turbulence: the Kraichnan spectrum in two dimensions and the Kolmogorov 
spectrum in three dimensions. 

The energy spectrum E{k) is a function that has the property that 



/ E{k)dk = {\u\^) . 
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(The convention is that < • • • > is normahzed space integral followed by long 
time average). The Kolmogorov spectrum for 3D turbulence is 

E{k)^CKi (e)^fc-t. 
The 2D Kraichnan spectrum is 

E{k)^CKr {v)^k-\ 

Here e = i/|Vup is the rate of dissipation of energy and rj = z^jVajp is the rate of 
dissipation of enstrophy. The spectra are supposed to be valid in a range of scales 
k € [ki, kd] where kd is the dissipation scale and is determined by viscosity and 
e (respectively viscosity and 77) alone. Their expressions arc then determined by 
dimensional analysis. We consider the Littlewood-Paley decomposition of the 
velocity associated to a moUifier (j>{£,). This is a smooth function in that 
is non-increasing, smooth, radially symmetric, satisfying (f){£^) — 1 for |^| < |, 

m - for lei > f . One sets V(o)(0 = " HO, ^(m)(e) = ^(o)(2-™0 and 

0(0 = / e-*(«-")$(z)dz, 

V'm(0= / e-^(«*)vl,(,„)(z)d^. 
The Littlewood-Paley decomposition is 

00 

where 

M(-oo) (x, t) = L-'^ $ (^1) u{x - y, t)dy, 

U{m){x,t) = L^'^ j ^ u{x ~ y,t)dy 

and L is a length (the integral scale). 

We define the Littlewood-Paley spectrum to be 

i?Lp(/c) = /c,-i(|ii(„)|2) 

for km-i < fc < km, m > 1 with km — 2™L~^. 

We start with the two dimensional Navier-Stokes equation 

{dt + u-V ~ iyA)uj ^ f 

with 

u{x,t)^—l j^uj{x-y,t)dy 
27r Jr2 \y\^ 

One can prove ([p4|) 
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Theorem 9 Assume that the source of vorticity in the Navier-Stokes equations 
has spectrum localized in the region of wave numbers fc < for some L > 0. 
Then there exists a constant C such that the Littlewood-Paley energy spectrum 
of solutions of two dimensional forced Navier-Stokes equations obeys the bound 



1 \ 6~ 

kd 



ELp{k) < Ck- 



for k > L^^. Here t^^ ||Vm||l=- >• 

The corresponding three-dimensional energy spectrum resuh requires a sig- 
nificant assumption: 

(IVmP) < oo. 



Denoting 



e^y{{\Vu\^)Y 

3 _1 

kd = i^^Heji 



and setting 



(^i'= J J \y'^i,){a)\\a\^\^io)ib)mdadb. 

we have (||): 

Theorem 10 Consider three-dimensional body forces that satisfy 

m = 

for all \k\ > and some C > 0. Consider solutions of the three dimensional 
Navier-Stokes equation that satisfy ? < oo. Then 

k \ 3 



ELp{k) <C^{t)^ k-^ ^ 
\kd 

holds for \k\> j^. 

Acknowledgment. Research partially supported by NSF DMS-9802611, by 
AIM, and by the ASCI Flash Center at the University of Chicago under DOE 
contract B341495. 



Bibliography 

1. U. Frisch, Turbulence: the legacy of A.N. Kolmogorov, Cambridge Univ. Press, 
(1995). 

2. K.R. Sreenivasan and R. Antonia, The phenomenology of small-scale turbu- 
lence Annu. Rev. Fluid Mech. 29 (1997), 435-472. 



12 



Navier-Stokes Equations and Fluid Turbulence 



3. A. Tsinober, Turbulence: beyond phenomenology, Lect. Notes. Phys. 511 
(1998), 85. 

4. C. Foias, C. Guillope, R. Temam, New a priori estimates for Navier-Stokes 
equations in dimension 3, Commun. Partial Diff. Eqns 6 (1981), 329-359. 

5. A. Majda, Vorticity and the mathematical theory of incompressible flow 
Comm. Pure Appl. Math. S39 (1986), 187-220. 

6. J.T. Bcale, T. Kato, A. Majda, Remarks on the breakdown of smooth solu- 
tions for the 3-D Euler equations, Commun. Math. Phys. 94 (1984), 61-66. 

7. P. Constantin, Geometric statistics in turbulence, SI AM Review 36 (1994), 
73-98. 

8. P. Constantin, A. Majda, E. Tabak, Formation of strong fronts in the 2D 
quasi-geostrophic thermal active scalar, Nonlinearity, 7 (1994), 1495-1533. 

9. P. Constantin, C. Fefferman, A. Majda, Geometric constraints on potentially 
singular solutions for the 3-D Euler equations, Commun. in PDE 21 (1996), 
559-571. 

10. D. Cordoba, On the geometry of solutions of the quasi-geostrophic active 
scalar and Euler equations, Proc. Natl. Acad. Sci. USA 94 (1997), 12769- 
12770. 

11. P. Constantin, An Eulerian-Lagrangian approach for incompressible fluids, 
preprint AIM (http://www.aimath.org/preprinte) 

12. J. Serrin, Mathematical principles of classical fluid mechanics, (S. Flugge, 
C. Truesdell Edtrs.) Handbuch der Physik, 8 (1959), 125-263, p.l69. 

13. M. E. Goldstein, Unsteady vortical and entropic distortion of potential flows 
round arbitrary obstacles, J. Fluid Meek. 89 (1978), 433-468. 

14. M. E. Goldstein, P. A. Durbin, The effect of finite turbulence spatial scale 
on the amplification of turbulence by a contracting stream, J. Fluid Meek. 98 
(1980), 473-508. 

15. J. C. R. Hunt, Vorticity and vortex dynamics in complex turbulent flows. 
Transactions of CSME, 11 (1987), 21-35. 

16. M. Brenner, P. Constantin, L. Kadanoff, A. Schenkel, S. Venkataramani, 
Diffusion, attraction and collapse, Nonlinearity 12 (1999), 1071-1098. 

17. P. Ronney, Some open issues in premixed turbulent combustion, in Modeling 
in Combustion Science^ J.Buckmuster and T.Takeno, eds.. Springer- Verlag, 
Berlin, 1995. 

18. P. Constantin, A. Kiselev, A. Oberman, L. Ryzhik, Bulk Burn- 



ing Rate in Passive- Reactive Diffusion, preprint math.AP/9907132 



http: / /xxx.lanl.gov/abs/math. AP /9907132 



19. L.N. Howard, Heat transport in turbulent convection, J. Fluid Mechanics 
17 (1964), 405-432. 

20. C. R. Doering, P. Constantin, Variational bounds on energy dissipation in 
incompressible flows HI. Convection, Phys. Rev E, 53 (1996), 5957-5981. 



Navier-Stokes Equations and Fluid Turbulence 



13 



21. B. Castaing, G. Gunaratnc, F. Hcslot, L. KadanofF, A. Libchabcr, S. 
Thomae, X.-Z. Wu, S. Zaleski and G. Zanetti, Scaling of hard thermal turbu- 
lence in Rayleigh-Benard convection, Journal of Fluid Mech. 204 (1989), 1-30. 

22. P. Constantin, C. Hallstrom, V. Putkaradze, Heat transport in rotating con- 
vection, Phys. D 125 (1999), 275-284. 

23. P. Constantin, C. Doering, Infinite Prandtl Number Convection, J. Stat. 
Phys. 94 (1999), 159-172. 

24. P. Constantin, The Littlewood-Paley spectrum in 2D turbulence, Theor. 
Comp. Fluid Dyn. 9 (1997), 183-189. 

25. P. Constantin, Q. Nie, S. Tanveer, Bounds for second order structure func- 
tions and energy spectrum in turbulence Phys. Fluids (in press) (1999). 



